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INNOVATIVE  METHODOLOGY  FOR  CANCELLING  CONTAMINATING  NOISE 
IN  TURBULENT  FLUID  FLOW  ENVIRONMENTS 


INTRODUCTION 

There  has  been  much  interest  in  the  last  few  years  concerning  the  origin  and  related  char¬ 
acteristics  of  the  wail  pressure  fluctuations  generated  at  the  wall  within  turbulent  flows.  For 
the  most  part,  such  investigations  have  been  concerned  with  the  coupling  phenomena  between 
the  solid  structure  (wall)  and  the  surrounding  fluid  which  may  result  in  undesired  structural 
vibration  and/or  radiated  acoustic  sound.  Such  fluid/structure  coupling  may  approach  reso¬ 
nance  conditions  wherein  fatique  of  the  solid  surface  material  may  occur  and  contaminating 
noise  radiated  to  the  farfield.  This  interaction  between  the  fluid  and  constraining  structure  pose 
specific  engineering  problems  to  the  Navy  in  terms  of  Sonar  and  towed  array  contamination 
and  the  detection  liability  resulting  from  radiated  acoustic  noise. 

Lately,  the  flow  induced  pressure  fluctuations  themselves,  have  proved  useful  in  furthering 
our  knowledge  of  the  fundamental  physics  of  turbulence.  Since  the  measurement  takes  place 
at  the  wall,  a  non-intrusive  examination  of  the  origins  of  turbulence  may  be  attempted  with 
no  disturbance  of  the  flow.  Such  basic  turbulent  characteristics  as  similarity  scaling,  length 
and  time  scales,  and  the  interaction  between  inner  and  outer  flow  regions  may  be  examined 
through  the  wall  pressure  field.  Recently,  Her(1986)  has  argued  that  the  Burst/Sweep  events  in 
the  near-wall  region  of  a  turbulent  flow  are  a  major  contributor  to  the  high  frequency  pressure 
fluctuations.  These  fluctuations  occur  within  the  inner  flow  region  of  a  boundary  layer  and  yet 
appear  to  scale  on  low  frequency  large  scale  motions  which  are  governed  by  the  outer  flow. 

Although  the  investigation  of  turbulent  flow  induced  pressure  fluctuations  is  highly  de¬ 
sirable,  obtaining  high  quality  measurements  has  proved  to  be  very  difficult.  The  reviews  by 
Willmarth(1975)  and  Leehey(1989)  have  discussed  in  some  detail  the  difficulties  associated 
with  resolving  the  frequency/ wavenumber  content  of  wall  pressure  fluctuations.  Both  reviews 
address  the  problems  associated  with  resolving  the  fine  scale,  high  wavenumber  portion  of 
the  spectrum  as  well  as  the  effects  of  signal  contamination  by  facility  generated  noise  at  low 
frequencies.  Accurate  resolution  of  this  low  frequency  end  of  the  wall  pressure  spectrum  has 
proved  elusive  in  almost  all  investigations  reported  in  the  literature  in  both  internal  and  ex¬ 
ternal  flow  configurations.  Willmarth  explains  that  the  sound  field  at  low  frequency,  either 
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flow  induced  or  facility  generated,  contaminates  low  frequency  data  and  is  most  difficult,  if  not 
impossible  to  eliminate  in  a  systematic  way  in  all  studies  done  to  date.  It  is  this  elimination 
of  low  frequency  noise  we  wish  to  address  in  this  work. 

BACKGROUND 

There  have  been  a  number  of  attempts  in  the  past  to  eliminate  or  reduce  unwanted  low 
frequency  noise  from  interferring  with  turbulent  flow  wall  pressure  fluctuation  measurements. 
The  early  methods  utilized  sound  isolators,  absorbers  and  polymeric  linings  to  damp  out  the 
contaminating  noise.  In  most  cases,  a  considerable  amount  (if  not  all)  of  the  low  frequency 
pressure  spectrum  was  eliminated  from  the  resulting  measurements.  Willmarth(1970)  intro¬ 
duced  a  new  technique  by  recognizing  that  typical  experimental  facilities  contain  local  source 
mechanisms  of  acoustic  sound  such  as  vanes,  diffusers,  bends  in  piping  and  recirculating  pumps. 
He  showed  that  a  finite  correlation  was  obtained  for  negative  time  delay  between  flush-mounted 
transducers  separated  in  the  streamwise  direction.  By  combining  this  temporal  quantity  with 
the  spatial  separation  between  two  sensors,  a  velocity  was  calculated  equal  to  the  sound  speed 
(relative  to  the  fluid  velocity)  in  the  fluid.  Willmarth  interpreted  this  quantity  as  facility 
generated  noise  and  subtracted  its  contribution  from  the  original  data. 

Wambsganss  and  Zaleski(1970)  introduced  a  signal  processing  scheme  which  was  later  re¬ 
fined  by  Wilson  et  al. (1979)  and  Horne  and  Hansen(1981).  This  method  utilized  three  pressure 
transducers  located  in  a  plane  perpendicular  to  the  flow  direction.  Using  simple  subtraction 
of  the  time-dependent  pressure  signals  and  then  cross-correlating  various  combinations  of  sub¬ 
tracted  signals,  a  means  of  separating  the  individual  outputs  due  to  turbulent  wall  pressure, 
background  acoustic  and  electronic  noise  could  be  resolved  respectively.  Although  innovative, 
these  temporal  subtraction  schemes  have  resulted  in  two  potential  sources  of  error.  The  first 
source  of  error  is  related  to  the  restriction  placed  on  the  contaminating  noise  which  is  assumed 
acoustic  in  origin  and  required  to  propagate  as  a  plane  wave.  This  requires  all  noise  to  originate 
in  the  far  field  and  doesn’t  provide  a  means  of  eliminating  higher  order  modes  of  propagation  or 
sources  of  noise  more  locally  generated.  The  second  source  of  error  is  related  to  the  fact  that  all 
previous  cancellation  schemes  involving  the  addition  and/or  subtraction  of  transducer  output 
signals  a  priori  assume  that  the  individual  transducer  sensitivities  are  equal  in  magnitude  and 
phase.  This  seems  trivial,  but  most  research  in  the  laboratory  deals  with  voltages  rather  than 
the  actual  pressures  themselves.  As  shown  by  Horne(1990),  a  finite  error  is  always  present  in 
typical  cancellation  schemes  related  to  the  ratios  of  the  sensitivities  of  the  transducers  used. 
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In  the  analysis  which  follows,  actual  pressure  signals  are  assumed.  Since  the  methodology 
involves  methods  of  statistical  processing  of  time  dependent  signals,  this  would  require  all 
sensor  signals  to  be  digitized  and  stored  on  tape.  Subsequent  application  of  the  cancellation 
technique  would  be  conducted  via  software  on  a  computer  after  all  appropriate  sensitivities 
have  been  applied.  Of  course,  if  the  individual  transducer  sentivities  can  be  made  equal  to 
each  other,  there  is  no  reason  the  technique  cannot  be  applied  in  the  laboratory  utilizing  a 
dual  channel  analyzer. 


CANCELLATION  METHODOLOGY 

The  techniques  described  previously  are  temporal  by  nature,  and  work  reasonably  well 
provided  the  sensor  sensitivites  are  equal  and  the  contaminating  noise  is  acoustic  in  origin  and 
propagates  as  a  plane  wave.  These  assumptions  are  restrictive  in  that  they  require  matched 
sensors  and  don’t  provide  a  means  of  cancelling  random  vibrational  energy  that  is  generated 
local  to  the  transducer  measuring  point.  The  method  to  be  presented  utilizes  the  concept 
of  a  “Wiener  filter”  (Wiener, 1949).  This  filter  was  designed  originally  to  cancel  wide  band 
frequency  signals  in  electrocardiography  and  the  broad-band  interference  in  the  sidelobes  of 
an  antenna  array.  In  its  original  definition,  the  Wiener  filter  assumes  that  a  signal-free  refer¬ 
ence  is  available  which  is  representative  of  only  the  contaminating  noise.  This  is  not  usually 
possible  in  most  turbulent  flow  situations.  However,  in  applying  the  filter  technology  to  the 
measurement  of  wall  pressure  fluctuations,  certain  assumptions  provide  a  unique  derivation  of 
a  set  of  correction  relationships  which  precisely  define  the  maximum  error  present  in  having  a 
contaminated  reference  signal  (a  non-signal  free  reference).  In  the  following  presentation,  the 
method  of  cancellation  is  formulated  together  with  an  innovative  error  analysis  applicable  to 
turbulent  channel  flow.  Afterwards,  some  representative  results  utillizing  the  methodology  will 
be  presented  and  discussed. 

The  measurement  scenario  can  be  demonstrated  by  the  schematic  illustrated  in  figure  1. 
As  shown,  there  are  two  contaminated  signals  represented  by  y\(t)  and  y2 ( ^ ) -  An  assumption 
is  made  that  the  contaminated  noise  is  present  in  each  signal  and  there  are  N  reference  signals 
available  for  cancellation  purposes.  No  a  priori  assumptions  are  made  about  the  sources  of 
the  noise  nor  how  it  propagates,  only  that  it  be  present  in  some  measurable  form  within  both 

signals.  Figure  1  shows  that  the  reference  signals  xi,i2,x3, . ,xyv  are  passed  through  two 

optimum  Wiener  filters  represented  by  hu,h2i, . ,hm  and  hu,h22, . ,hN2.  These  filtered 


3 


signals  awe  then  subtracted  from  the  original  signals  y\(t)  and  y2(t)  to  obtain  the  desired  “noise- 
free”  signals,  Z\(t)  and  Z2{t).  The  goal  here  is  to  minimize  the  mean  square  output  and  obtain 
the  auto  and  cross-spectral  properties  of  z\(t)  and  22(0  respectively. 

In  the  following,  all  calculations  are  performed  on  Fourier  transformed  quantities  via  the 
following  relationship. 

X*  =  Xk(f)  =  f"  xk(t)e-^^dt  (1) 

The  variable  xk(t)  is  real,  but  all  Fourier  transformed  quantities  of  xk(t)  are  assumed 
complex  and  implicit  functions  of  frequency,  /.  A  superscript  ‘  *’  denotes  a  complex  conjugate 

of  a  Fourier  transform,  hnd  the  brackets  (braces)  < . >  represent  the  statistical  expectation 

of  all  ensembles  evaluated.  Bearing  this  in  mind,  the  following  (boldface)  vectors  can  be  defined. 


X  =  [Xi  X2  X3  ...  XN]  (2a) 

Y  =  [Yx  Y2]  (2b) 

1  =  [Zx  Z2)  (2c) 


'  fin 

h\2  ‘ 

h2\ 

h22 

H  = 

/131 

h 32 

■  hN\ 

hp/2 . 

(2d) 


The  block  diagram  indicated  in  figure  1  can  be  represented  by  the  following  vector  equation, 
using  equation  set  (2). 


Z  =  Y  -  X  H  (3) 

We  can  make  the  following  matrix  definitions  for  the  auto  and  cross-spectral  properties  for  the 
signals,  X),  Y,,  and  Z{.  (Note,  that  the  cross  spectrum  of  the  Fourier  transforms  of  x<(t)  and 
Xj(t)  is  defined  as,  GXi X)  =  <  Xm  X  >.) 
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’  Gi,x, 

GXix3 

•  •  •  Gx\XH 

Gxjxx 

Gx3X3 

•  •  •  Gx3XN 

=  <  X*  X  >  = 

Gx3X\ 

Gx3X3 

•  •  •  Gx3xn 

■GXsx  1 

GxnXj 

■  •  -  GXnxn  • 

R  =  <  Y*  Y  > 


^VlVl  ^»l»J 

Gy3yx  ^V2»a  . 


(4a) 


(46) 


P  =  <  Zm  Z  > 


^2l*J 

.^2321  ^222a 


(4c) 


The  auto/cross-spectral  properties  of  the  output  signals  z\  and  z2  are  desired  for  which  the 
following  manipulation  may  be  made  with  equation  (3). 

P  =  <  Z*Z  > 


=  <(Y  -  XH]*  [Y  -  XH]  > 

=  <  [Y"  -  H*X*]  [Y  -  XH]  >  (5) 

Note,  the  definition  of  a  Hermitian  conjugate,  ( A{jBjk)m  =  B^A^  has  been  used  in  the 

above  derivation.  Performing  some  algebra  with  equation  (5),  results  in  the  following, 

P  =  R  -  [<  Y*X  >  H]  -  [H*  <  X*Y  >]  +  [H'QH]  .  (6) 

The  following  matrix  may  be  defined  which  corresponds  to  the  cross-spectral  properties  of  the 

Xi  and  Yi  signals. 


V  =  <  Y“  X  > 


Gyixi  GViXj 
.Gyan  GVlXj 


GyixN 

Gy7xN  . 


(T) 


Note,  that  V*  =  [<  (Y“X)*  >]  =  (<  X*(Y*)*  >]  =  <  X*Y  >.  Substituting  V  and  V* 
appropriately  into  equation  (6)  results  in  the  following  relation. 

P  =  R  -  VH  -  H*V*  +  H*QH  (8) 
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By  substituting  an  expression  for  the  optimum  Wiener  filter,  P  can  be  expressed  appropriately 
in  terms  of  the  measured  signal  quantities,  R,  V,  and  Q.  It  can  be  shown  that  the  optimum 
Wiener  filter  which  minimizes  the  mean  square  output  power  of  equation  (8)  is  given  by, 

H  =  (Q-1)*  V*  .  (9) 

Substituting  equation  (9)  and  its  conjugate  into  equation  (8)  results  in, 

P=<Z*Z>=R  -  VQ_1V*  .  (10) 

Equation  (10)  represents  the  final  result.  Referring  back  to  the  definitions  of  R,  V  and  Q, 
one  may  note  that  these  three  vectors  are  in  terms  of  auto  and  cross-spectral  properties  of 
given  measurable  signals.  Also  note,  that  equation  ( 10)  eliminates  the  noise  identically  if  the 
reference  signals  Xj(t)  are  truly  "signal  free”.  That  is,  the  reference  signals  contain  only  the 
noise  which  is  present  in  the  original  signals.  If  the  reference  signals  are  not  truly  "signal  free”, 
then  equation  (10)  will  only  give  the  "best”  estimate  possible.  Under  certain  situations,  exact 
estimates  of  the  inherrent  error  may  be  made,  as  will  be  shown  later. 

SPECIAL  CASE  OF  ONE  REFERENCE  SIGNAL 

Suppose  there  is  only  one  reference  signal.  Then  Q  reduces  to  a  scalar  defined  by, 

Q  =  <  X*X  >  =  GXlXl  =  Gzx  ,  (11) 

and 

Q~l  =  1  /Gxx.  (12) 

Substituting  equation  (12)  into  equation  (10),  and  recognizing  that  V  and  V*  are  two-element 
vectors,  the  following  relations  are  obtained  for  the  auto  and  cross-spectrums  of  the  signals 
zi{t)  and  Z2(0- 

GZln  —  Gv jyj  —  -p;  GylX  GV}X  (13a) 

=  GVl  y,  —  p;  GylX  GyoX  (136) 

XT 

G  Z7*\  =  Gy7  y,  —  Gy^  X  G  y,  r  (13c) 

^  XX 

G Z7Z7  =  Gy7y2  —  —  Gy7l  G ylX  (13d) 

O  XX 


Note,  the  order  of  the  indices  in  equation  set  (13)  is  imperative  and  equations  (13b)  and  (13c) 
are  for  all  intents  and  purposes,  equivalent.  If  the  signals  jq  and  yi  axe  of  the  following  form, 


-  31  +  *  , 

(14a) 

Vi  =  -53  +  *  > 

(146) 

and  Si  and  32  are  neither  correlated  with  each  other  or  with  x ,  the  reader  may  verify  for  himself 
that  equation  (10)  eliminates  the  noise  exactly  given  that  the  desired  signals  are  GSI j, ,  Gtlt3 
and  respectively. 

SINGLE-POINT  APPLICATION 

Figure  2  shows  a  typical  transducer  alignment  for  single  and  two- point  cross-spectral  mea¬ 
surements  of  the  pressure  field  at  the  wall  within  a  boundary  layer  flow  such  as  a  two  dimen¬ 
sional  channel.  No  restriction  is  attributed  to  the  noise  signal  as  to  its  frequency  content  or 
propagational  characteristics.  The  total  signal  from  each  transducer  is  expected  to  be  of  the 
form  indicated  in  equation  (14)  for  which  s(f)  represents  the  turbulence  and  x(t)  represents 
the  noise.  Based  on  this  composition,  the  following  assumptions  are  adherred  to. 

(a)  The  turbulence  at  any  location  is  not  correlated  with  the  noise  signal  at  its 
own  location,  nor  with  that  at  any  other  location. 

(b)  The  noise  signal  at  “every”  location  is  either  fully  correlated  or  has  finite 
measurable  correlation  with  that  at  any  of  the  other  locations. 

(c)  It  is  postulated  that  the  ratio  l/Ax  >»>  1  such  that  the  turbulence  at 
location  ‘1’  is  uncorrelated  with  the  turbulence  at  any  of  the  other  locations. 
The  same  is  not  true  between  transducer  2/3,  2/4,  and  so  forth. 

From  the  definition  of  a  “Wiener"  filter  and  having  one  reference  signal  for  cancellation 
purposes,  equation  set  (13)  can  be  recast  in  terms  of  i  and  j  as, 

GstZi  ~  GyiVi  —  ~  GyiX  Gy.x  ,  (15a) 

GziZ,  =  GylVj  —  ~  GVix  Gy  J.  ,  (156) 

U'  XX 
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where  i  =  2  and  j  =  3,  4,  5,  N.  Suppose  we  make  the  following  definitions  for  the  signals 
for  each  transducer. 

Vi{t)  =  x(t)  (16a) 

2/2  (0  =  x(t)  +  s2(t)  (166) 

yj(t)  =  x{t)  +  Sj(t),  for  j  =  3,4 . ,N  (16c) 


In  other  words,  yi(t)  will  be  used  as  a  reference  signal.  The  quantities  desired  are  the 
spectral  functions,  G,2B2,  and  Gi7Bj  for  j  =  3,  4,  5,  6,  ....  ,  N.  Note,  nothing  has  been  said 
about  the  origins  of  the  noise  signal,  only  that  it  be  the  ‘same’  at  each  location.  Recognize  that 
a  specific  type  of  noise  represented  by  acoustic  noise  that  propagates  as  a  plane  wave  satisfies 
the  aforementioned  assumptions  identically  in  a  cross-sectional  plane  of  a  channel.  In  the  flow 
direction,  this  is  not  necessarily  true.  However,  if  the  acoustic  noise  originates  considerably 
in  the  farfield,  for  all  intents  and  purposes,  the  noise  at  each  transducer  location  is  identical. 
Using  equation  (15a)  we  have  the  following. 


G 


'yiyi 


+  GIX]  -  — —  G2XX 

G  sjtj 


(17) 


Similar  results  are  obtained  for  the  cross-spectral  density  function,  GZiZ/ .  Thus,  for  a  signal 
free  reference,  the  noise  can  be  eliminated  exactly.  This  brings  us  to  an  interesting  point.  In 
almost  any  experimental  facility,  it’s  virtually  impossible  to  obtain  a  signal  representative  of 
the  noise  only.  If  transducer  T’  is  just  another  wall  pressure  transducer,  it’s  not  signal  free; 
hence  we  have  a  contaminated  cancellation  reference.  Assuming  that  the  signal  y\(t)  can  be 
represented  by; 

yi(t)  =  jj(t)  +  x(t)  , 


then  equation  (15a)  takes  the  following  form  with  Gxx  replaced  with  GVlVl 


-  [G33tj  +  Gxx]  - 


Glr 


G,iBj  +  Gz 


(18) 


At  this  point,  an  interesting  observation  may  be  made.  For  a  fully  developed  turbulent 
channel  flow,  it  may  be  assumed  that  all  turbulence  and  related  statistical  quantities  are  inde¬ 
pendent  of  the  spatial  dimension  in  the  flow  direction.  In  addition,  this  type  turbulent  flow  has 
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“homogeneous”  statistical  properties,  particularly  as  represented  by  auto-spectral  densities. 
Hence,  G ,,3,  =  G3333  =  G3i3)  for  all  *j\  Therefore,  in  equation  (18)  can  be  replaced 

with  G3 3,3.  Letting  Q2  represent  the  noise  to  signal  ratio  Gxx/G3332,  and  performing  some 
algebra  results  in  the  following  equation. 

=  g^[Wa°~]  <19) 

To  summarize  briefly,  equation  set  (15)  has  been  shown  to  correctly  resolve  the  desired 
turbulent  signal  for  transducer  ‘2’  for  a  ‘signal-free’  reference  at  location  ‘1’.  However,  if 
transducer  l’s  signal  also  contains  turbulence,  equation  (19)  results.  Three  cases  can  be  defined 
for  values  of  £*2  • 

case(l)  a2  <<<<<  1  — ♦  0  (Gxx  =  0) 

Gz2z2  —  G 3333  Thus,  if  no  noise  is  present  in  sensor  2,  the  use  of 
the  cancellation  methodology  introduces  no  additional  error. 

case  (2)  q2  =  1  (Gxx  =  G3232) 

Gt3Z3  =  |  G„3 j2,  or  a  signal  approximately  1.76  dB  too  high. 


case  (3)  £*2  >>>>>  1  — *•  00  (noise  very  large) 

GZJZJ  =  2  G3332,  or  a  signal  that  at  most,  is  only  3  dB  larger  than 
it  should  be. 


Thus  we  see  from  these  3  cases,  equation  (19)  either  introduces  no  additional  contami¬ 
nation,  or  that  the  desired  result  is  obtained  with  at  most,  an  error  of  3  dB.  Knowledge  of 
the  level  of  noise  is  not  know  a  priori ,  but  some  insight  may  be  obtained  by  looking  at  the 
coherence  between  transducers  1  and  2.  The  coherence  is  defined  by, 

_  I  G.»,  I* 

7“*'  aw  g„„  ’ 


or, 


g| r 

(G3332  +  G«)2  ’ 
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2 

*3*1 


(  <*1  )2 
(1  +  q2  )2 


(20) 
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The  coherence  ranges  from  0  to  1  with  7^  =  0,  corresponding  to  case  (1)  and  7*  =  1, 

corresponding  to  case  (3).  Equations  (18)  and  (19)  comprise  a  new  process  for  cancelling  noise 
from  turbulent  wall  pressure  fluctuation  measurements.  The  first  applies  an  innovative  method 
for  cancelling  unwanted  noise  using  a  reference  signal,  and  the  second  attempts  to  correct  for 
this  reference  signal  being  contaminated,  assuming  the  noise  is  the  same  at  each  location. 

So  far,  we’ve  only  investigated  the  auto-spectral  signals.  The  question  can  well  be  asked 
whether  the  methodolgy  (and  subsequent  correction)  can  be  similarly  defined  for  the  two-point 

spectral  statistics  obtained  from  equation  (15b),  such  as  G,2Jf  with  j  =  3,  4,  . ,  .V.  The 

answer  is  yes,  with  an  interesting  result. 

MULTI-POINT  APPLICATION 

In  order  to  address  the  issue  of  cross-spectral  properties  we  examine  equation  (15b)  with 
i  =  2  and  j  =  3  as  an  example.  This  equation  takes  the  form, 

G2g*3  =  Gy2y 3  ~~  p;  Gy3yi  ,  (21) 

where  we  have  substituted  y±(t )  =  3i(t)  4-  x(t)  since  transducer  1  is  assumed  to  be  a 

contaminated  reference  signal.  We  may  note  that  according  to  the  original  assumptions,  the 
following  equalities  may  be  made. 


G  yjy3 

=  G,^,3 

+ 

Gxx 

(22  a) 

G  yjy! 

=  GXI 

(22  b) 

Cm 

=  Gxx 

(22c) 

G  yiyi 

=  G4lJ, 

+ 

Gxx 

(22d) 

Substituting  these  relations  into  equation  (21)  results  in  the  following. 

Gz3z3  =  G,  +  Gs,,,  —  T  (23) 

(1  +  Qi) 

Note  that  the  turbulent  cross-spectrum  G,2,3,  is  complex  whereas,  the  auto  spectrum  G,,,,  is 
real.  Hence,  the  cancellation  method  incurrs  an  error  in  the  “real”  part  of  the  cross-spectrum 
only.  This  is  not  unexpected  since  we  have  assumed  the  noise  to  not  comprise  a  convective 
component  between  transducers.  Again,  Q\  is  defined  as  the  ratio  of  the  noise  to  the  level 
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of  turbulence  (GXx/GSltl ).  It  can  be  shown  for  the  assumptions  outlined  above,  that  the 
correction  equations  (19)  and  (23)  would  take  the  following  form  for  ”k”  reference  transducers. 


G 


*2*2 


r  1  +  (k  4-  1)  oi2 

°l3J3l  1  +  ka2 


(24) 


'  *2  *3 


=  G 


*2*3 


+  G 


>i*i 


( <*i  ) 


(l  +  k  cm) 


(25) 


Equations  (24)  and  (25)  show  that  as  the  number  of  reference  transducers  available  for  cancel¬ 
lation  increases,  very  little  correction  is  needed  for  the  signals  being  contaminated,  as  would 
be  expected  intuitively. 


Note,  we  can’t  emphasize  enough  that  this  correction  as  stipulated  by  equations  (23) 
through  (25),  are  accurate  only  to  the  extent  that  the  noise  contained  in  all  signals  is  identical. 
This  a  priori  requires  the  noise  to  be  generated  in  the  farfield  at  very  large  wavelengths, 
neglecting  a  convective  component  in  the  streamwise  flow  direction.  For  locally  generated  noise 
sources,  or  those  which  propagate  at  other  modes  than  that  of  a  plane  wave,  more  reference 
signals  would  be  required  in  order  to  correct  for  the  convective  properties  of  the  noise  source 
term. 


EXPERIMENTAL  ARRANGEMENT 

In  order  to  test  the  methodology  presented  above,  a  plane  Poisuelle  flow  was  investigated  in 
an  attempt  to  obtain  high  quality  measurements  of  the  wall  pressure  field.  These  measurements 
were  made  in  a  457  mm  by  25  mm  (18:1  aspect  ratio)  rectangular  channel  flow  facility  located 
at  the  Naval  Research  Laboratory  with  water  as  the  working  fluid.  The  channel  operates  in  a 
blowdown  mode,  wherein  pressurized  air  is  added  to  an  upstream  reservoir  of  water  at  such  a 
rate  as  to  maintain  an  overall  driving  pressure  throughout  a  particular  experimental  run.  The 
test  section  contains  a  number  of  large  plates  within  which  smaller  test  windows  or  transducer 
plates  may  be  arranged  in  many  different  configurations.  A  schematic  of  this  facility  is  shown  in 
figure  3.  The  measurement  location  was  4.1m  downstream  from  the  beginning  of  the  developing 
turbulent  flow.  The  ratio  of  the  downstream  distance  of  the  measurement  point  to  the  half¬ 
height  of  the  channel  was  320.  This  value  was  within  the  range  for  which  fully-developed 
turbulent  channel  flow  can  be  expected,  as  shown  by  Hussain(  1975). 

To  validate  the  existence  of  a  generic  turbulent  flow  within  the  channel,  representative 
measurements  of  the  streamwise  velocity  were  obtained  for  a  range  of  Reynolds  numbers. 
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=  Uh/v,  and  are  shown  in  figure  4.  (In  this  experiment,  U  represents  the  centerline  velocity 
in  the  channel.)  This  data  is  plotted  in  log-law  coordinates  and  shows  excellent  agreement  with 
the  classical  equation  derived  by  Spalding(1961)  as  represented  by  the  solid  line. 

Five  pressure  transducers  (ENDEVCO  Model  8514-10)  were  flush  mounted  within  the 
channel’s  test  section  in  both  spanwise  and  streamwise  directions  as  depicted  in  figure  2.  The 
nominal  sensitivity  of  the  transducers  was  approximately  -227  dB  (re  1.0  ^f)  and  the  fre¬ 
quency  response  was  flat  out  to  140A: Hz.  The  active  area  of  the  transducer  face  is  equal  to 
0.5mm  (Galib  and  Zandina,1984)  giving  a  channel  half- height  to  sensor  diameter  ratio,  h/d,  of 
25:1  and  a  viscous  scale,  d+  =  duT/v,  in  the  range  20  to  40.  The  signals  were  low-pass  filtered  to 
prevent  aliasing,  then  digitized  and  stored  on  magnetic  tape  for  post-processing.  Software  was 
written  to  apply  appropriate  transducer  sensitivities,  calculate  Fourier  transforms  and  apply 
the  cancellation  methodology  presented  above. 

The  data  set  available  for  the  present  analysis  was  limited  as  to  the  total  number  of 
points  in  time  due  to  a  finite  volume  of  water  present  for  each  experimental  run.  According  to 
Bendat  and  Piersol(1980),  the  accuracy  in  measuring  spectral  statistical  functions  is  inversely 
proportional  to  Np,  the  number  of  ensemble  averages  obtained.  Higher  resolution  in  spectral 
space  requires  a  larger  transform  size  (larger  set  of  time  points),  resulting  in  an  overall  lower 
number  of  ensemble  averages  for  a  fixed  time  series.  The  spectra  presented  here,  were  obtained 
as  a  composite  of  two  Fourier  analyses;  one  with  a  resolution  of  1.95 Hz  and  one  at  a  resolution 
of  0.244 Hz,  with  approximate  spectral  accuracies  of  3.5%  and  9.5%  respectively. 

RESULTS  AND  DISCUSSION 

The  distance  /  separating  transducers  1  and  2  was  chosen  to  be  four  times  the  channel's 
half-height  h  such  that  the  turbulence  from  transducers  2  thru  5  had  zero  correlation  with  that 
obtained  from  the  reference  transducer  (transducer  1).  The  only  signal  components  remaining 
coherent  were  those  due  to  acoustic  and/or  facility  generated  noise.  Figure  5  presents  the 
power  spectral  density  obtained  from  transducer  2  for  the  Reynolds  number  R h  =  25,000. 
In  this  figure,  results  are  presented  utillizing  both  equation  set  (15)  as  well  as  the  further 
correction  obtained  from  equations  (19)  and  (20)  due  to  a  contaminated  reference  signal.  The 
un-corrected  total  sensor  output  indicates  that  the  channel  has  low  frequency  noise  which  peaks 
at  frequencies  of  1,  5  and  10  Hz  respectively.  This  noise  is  believed  to  be  caused  by  vortical 
separations  in  the  upstream  turbulence  management  section  and/or  posssible  contributions 
from  a  low  frequency  bulk  pressure  wave  resulting  from  the  blow-down  operational  mode  of 
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the  channel.  In  either  case,  the  results  indicate  that  more  than  40 dB  cancellation  is  obtained 
in  this  contaminated  region  with  excellent  resolution  at  higher  frequencies  where  limited  or  no 
contamination  exists.  The  line  representing  partial  correction  is  a  result  of  applying  equation 
set  (15)  without  the  correction  from  equation  (19).  Note  the  solid  line  representing  the  fully 
corrected  spectra  merges  smoothly  with  the  rest  of  the  spectrum  but  the  partially  corrected 
line  has  a  disjointed  connection  at  approximately  10  Hz.  Also  note  the  slight  waviness  of  these 
two  corrected  lines.  This  is  a  result  of  the  lower  spectral  accuracy  obtained  by  combining  two 
Fourier  analyses.  In  order  to  resolve  these  low  frequencies,  a  much  larger  transform  size  was 
required  resulting  in  a  lower  ensemble  average. 

Figure  6  shows  a  representative  plot  of  the  coherence  7 2j4i  and  figure  7  shows  the  max¬ 
imum  error  present  over  the  frequency  range  investigated.  This  error  can  be  represented  by 
rearranging  equation  (19)  as, 

$ 

error  =  10 log  =  IO/05 

As  figure  7  indicates,  the  maximum  error  present  equals  3 dB  which  may  also  be  interpreted 
from  figure  5  by  the  constant  3 dB  separation  between  the  lines  representing  the  partially  and 
fully  corrected  results  at  low  frequencies.  It  should  be  noted  that  in  most  cases  within  the 
laboratory,  calibrating  transducers  can  at  most  obtain  an  accuracy  of  approximately  0.5  to  1.0 
dB  corresponding  to  measuring  a  minimum  coherence  between  0.02  and  0.06.  This  is  at  the 
extreme  limitation  of  the  capabilities  of  any  dual  channel  analyzer,  hence  transducer  calibration 
appears  to  govern  the  limitations  on  applying  the  corrections,  rather  than  the  methodology 
itself. 

Figure  8  shows  the  representative  results  for  the  single  point  wall  pressure  spectra  at  three 
representative  Reynolds  numbers.  Chase(1980)  and  Panton(1980)  have  proposed  that  the  low 
frequency  part  of  the  wall  pressure  spectrum  should  drop  off  proportional  to  the  square  of 
frequency.  This  implies  a  d.op  of  20  dB  (or  more)  over  a  decade  in  frequency.  Figure  8 
indicates  an  approximate  /2  dependence  over  the  frequency  range  0.1  Hz  to  2.0 Hz.  Although 
this  trend  at  low  frequency  appears  to  be  validated,  this  result  is  of  a  preliminary  nature  and 
more  data  is  required  at  lower  frequencies  to  establish  this  relationship  more  definitively. 

As  an  illustration  of  the  necessity  of  using  cancellation  in  wall  pressure  experiments,  figure 
9  shows  a  typical  result  of  the  two-point  space-time  correlation  function  obtained  from  trans¬ 
ducers  2  thru  5  for  one  Reynolds  number.  One  may  note  that  the  three  solid  lines  indicated 


1  +  2a 


1  +  a 


(26) 
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at  the  top  of  the  figure  correspond  to  the  direct  correlations  with  no  cancellation  or  subse¬ 
quent  correction.  As  shown,  these  results  are  typical  of  wail  pressure  measurements  which  are 
contaminated  at  low  frequencies  by  high  intensity  noise  sources  of  extremely  long  wavelength. 
As  figure  9  indicates,  and  as  first  shown  in  figure  5,  the  facility  generated  noise  below  10  Hz 
completely  swamps  the  correlation  function,  giving  no  useful  information.  On  the  other  hand, 
by  applying  cancellation  together  with  the  additional  correction  for  a  non-signal  free  reference, 
the  more  definitive  results  shown  in  the  bottom  two-thirds  of  figure  9  may  be  resolved.  The 
successive  humps,  or  maximum  correlations  which  decay  in  magnitude  as  they  are  displaced  in 
time,  are  typical  of  a  convective  pattern  of  turbulence  which  rapidly  loses  coherence  as  may  be 
represented  by  a  turbulent  channel  flow. 

CONCLUSIONS 

A  method  has  been  presented  for  the  cancellation  of  low  frequency  acoustic  modes  which 
may  be  masking  significant  information  about  turbulent  wall  pressure  fluctuations.  The  method 
makes  use  of  certain  simplifying  assumptions.  Among  these  is  the  assumption  that  the  coher¬ 
ence  length  of  the  contaminating  mode  is  much  larger  than  that  of  the  turbulent  wall  pressure 
fluctuations  themselves.  The  method  also  incorporates  the  assumption  that  the  turbulence 
is  homogeneous  in  horizontal  planes  and  is  therefore  restricted  to  application  within  bound¬ 
ary  layer,  pipe  and  channel  flows.  These  assumptions  lead  to  a  two  part  correction  scheme 
in  which  first,  a  corrected  spectrum  is  found  using  a  standard  least  mean  squares  algorithm. 
Afterwards,  utilizing  the  assumption  of  homogeneity,  a  correction  factor  computed  from  the 
measured  coherence  is  applied  resulting  in  the  fully  corrected  auto  and  cross-spectral  density. 
It  is  also  shown  that  there  can  be  at  most  only  a  3  dB  difference  between  the  partially  corrected 
spectrum  and  the  fully  corrected  one.  This  3  dB  correction  is  obtained  in  the  limit  of  a  very 
large  noise  to  signal  ratio. 

Experiments  were  conducted  to  demonstrate  the  utility  of  the  method  in  the  case  of  fully 
developed  turbulent  channel  flow.  Wall  pressure  results  at  these  low  frequencies  for  internal 
flow  geometries  is  virtually  non-existent  within  the  literature.  This  new  cancellation  scheme 
may  provide  new  data  which  may  help  validate  current  theories  about  the  very  low  frequency 
behavior  of  the  wall  pressure  spectrum.  In  addition,  the  capability  of  resolving  two-point 
spectral  statistics  at  these  low  frequencies  is  of  particular  interest  to  the  Navy.  For  it  is  within 
this  low  frequency  range,  that  towed  array  topology  and  sonar  platforms  incur  maximum 
contamination  from  the  inherent  turbulent  flow  they  are  typically  subjected  to. 
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Fig.  1  —  Schematic  of  cancellation  method 
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Fig.  2  —  Transducer  arrangement  for  single  and  multi-point  spectral  statistics 
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Fig.  4  —  Universal  velocity  profile  for  various  Reynolds  numbers 
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Fig.  5  —  Power  spectral  density  at  Rh  =  25,000 
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Fig.  6  —  Coherence  between  the  two  spanwise  transducers 
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Fig.  8  —  Fully  corrected  power  spectral  density  for  various  Reynolds  numbers 
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